Partly noiseless encoding of quantum information in quantum dot arrays against 

phonon-induced pure dephasing 



N£> 
O 
O 



A. GrodeckfEl and P. Machnikowski 
Institute of Physics, Wroclaw University of Technology, 50-370 Wroclaw, Poland 

We show that pure dephasing of a quantum dot charge (excitonic) qubit may be reduced for 
sufficiently slow gating by collectively encoding quantum information in an array of quantum dots. 
We study the role of the size and structure of the array and of the exciton lifetime for the resulting 
total error of a single-qubit operation. 
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I. INTRODUCTION 

The idea of implementing quantum information pro- 
cessing schemes on charge states in quantum dot (QD) 
systemsi^ has recently been supported by an exper- 
imental demonstration^. On the other hand, it has 
been proposed that quantum bits encoded in electron 
spins (which have much longer dephasing times) may 
be optically controlled by inducing charge dynamics 
within various schemes exploiting selection rules and 
Pauli blockin g 4 ' 5 ' 6 ' 7 . Moreover, the demonstration of a 
QD-based photodiode operating in the quantum coherent 
regime^ shows great promise for quantum optoelectronic 
applications. These experimental achievements and the- 
oretical proposals have driven a lot of investigation on 
coherence properties of QD systems, in particular under 
external driving of carrier dynamics. 

Limitations to the coherent control of confined carriers 
in QDs result, on one side, from finite lifetime of charge 
states used for information encoding; one restriction is 
imposed, for instance, by radiative decay (on nanosec- 
ond time scale?) of confined cxcitons that seemed to be 
the most natural charge qubit systems due to relative 
ease of ultrafast optical driving 2 * On the other side, 
optical experiments^ have demonstrated that pure de- 
phasing processes destroy coherence of the system state 
to a large extent within a few picoseconds after creating 
the state. This has been attributed to carrier-phonon 
interactio n 11 ' 12 and interpreted in terms of spontaneous 
relaxation of the lattice after an optically induced change 
of the confined charge distribution, accompanied by emis- 
sion of phonon wave packetsi&ii, which creates a kind of 
which way trace in the crystal^ 5 - and results in dephasing. 
The same phonon response to the charge evolution leads 
also to dephasing in optical spin control scheme d 5 ' 16 . 

The phonon-induced dephasing effect may be avoided 
by driving the system evolution adiabatically with re- 
spect to phonon dynamics 5 ^ 7 ", which requires long driv- 
ing times and leads to growing probability of exciton de- 
cay during the operation. Optimal driving conditions 
correspond to a tradeoff between these two dephasing 
effects 17 . As a result, maximum fidelity of a coherent 
quantum control operation on a given system is limited. 

One way to partly circumvent this limitation is op- 
timizing the pulse shap o 18 ' 19 . Another solution might 
be to modify the properties of the phonon reservoir 
by placing the dot near the crystal surface 2 ^ or in an 



acoustic cavity 21 . An alternative approach is to over- 
come the physical restrictions imposed by dephasing pro- 
cesses on the quantum-logical level. If the reservoir 
acts collectively on a register composed of a number of 
charge qubits then subspaces exist which are safe against 
decoherence 2 ^ 2 ^ 4 * 2 ^. In principle, this idea can be used 
to implement robust logical qubits on QD array a 1 ' 26 . In 
practice, however, protection against decoherence based 
on collective interaction with the environment (closely re- 
lated to the subradiance effect) may be difficult to achieve 
in these artificial systems since it requires that the tran- 
sition energies in all QDs are the same. Moreover, the 
phonon-induced pure dephasing effect involves a range of 
modes selected by the system geometry in such a way 
that reservoir correlation length roughly corresponds to 
the single QD size which apparently precludes any col- 
lective interaction with a number of QDs. 

In this paper we show that, in spite of these difficul- 
ties, partly noiseless encoding against the pure dephasing 
effect is possible. We focus our discussion on the case 
of an excitonic charge qubit. In analogy to the original 
idea of noiseless encoding 2 ^, we will identify two states 
of an array of quantum dots which span a subspace (log- 
ical qubit) that, under appropriate driving conditions, 
allows performing single-qubit operations at reduced er- 
ror, i.e., with relatively little lattice response, compared 
to switching of a single excitonic physical qubit. Physi- 
cally, the collective encoding approach in this case relies 
on the fact that the pure dephasing effect is of dynami- 
cal character and results exclusively from nonadiabatic- 
ity of the drivingi 7 . For a given speed of rotation in 
the qubit space, it involves only slow modes that cannot 
follow the evolution of charge distribution, while faster 
modes follow adiabatically and reversibly. In this way, 
only a long wave length part of all the coupled modes is 
relevant for the dynamics. Since the range of these slow 
modes is selected dynamically, independently of the sys- 
tem geometry, it may correspond (under suitable driving 
conditions) to wave lengths larger than the register (QD 
array) size, thus allowing for collective interaction. 

In the present proposal, the logical qubit is defined on 
a system of N — 2™ quantum dots (physical qubits) in 
such a way that the logical |1) state corresponds to the 
presence of an exciton in each of the appropriately cho- 
sen N/2 dots, while in the logical |0) state the other half 
of the dots are occupied. With the correct definition of 
the logical qubit, the excitations of the phonon reservoir 
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induced by the evolution of charges in all QDs interfere 
destructively in the long wave length sector, leading to a 
decrease of dephasing. The proposed scheme works only 
for relatively slow evolution and the resulting fidelity gain 
is again restricted by finite lifetimes. We will study how 
much the overall error may be reduced depending on the 
exciton lifetime and the strength of other decoherence 
processes. It should be stressed that the scheme pro- 
posed here does not assume that the individual dots in 
the register are identical. 

The paper is organized as follows. In Sec. [H] we intro- 
duce the model of the QD array implementing the logical 
qubit. Next, in Sec. IHII we present the general method 
for describing the effects of phonon perturbation on an 
arbitrary operation on a logical qubit. Sec. IIVI contains 
the results, discussing the error reduction for various lay- 
outs of the logical qubit. Sec.lVlconcludes the paper with 
final remarks. 



II. MODEL 

As a model of the physical QD array implementing the 
logical qubit, we consider N quantum dots, each approxi- 
mated by a two-level system corresponding to the absence 
or presence of an exciton with a fixed polarization, with 
transition energies e;, I = 1, . . .N. The wave functions of 
carriers confined in different dots do not overlap so that 
no phonon-assisted transfer is possible. We will be inter- 
ested in driving the system with long, i.e., spectrally nar- 
row pulses so that excitations of optical phonons may be 
excluded. If the electron and hole wave functions in each 
dot overlap the piezoelectric coupling is very weak due 
to charge cancellation- 11 . Therefore, we assume that the 
confined carriers interact only with longitudinal acoustic 
phonons via deformation potential coupling. The Hamil- 
tonian of this semiconductor system is then 

H sc = ^ e„a^a„ + ^ ^ fc &j> fc (1) 

n k 
+ a l a nfn(k) (b{ + 6_fc) , 

n,fe 

where a\ , a n are fermionic creation and annihilation op- 
erators for the exciton in the nth dot, b k , bk are bosonic 
operators for phonon modes, LOk = ck are the correspond- 
ing frequencies (c is the speed of sound), and /„(fe) are 
the carrier-phonon coupling constants for an exciton in 
the nth dot with the symmetry f n (k) — /*(— k). As- 
suming, for simplicity, that the exciton state in the nth 
dot is described by a product of electron and hole wave 
functions ipii (i — r n ), where r n is the position of the 
dot, the coupling constants may be written as^ 
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Here p c denotes the crystal density, v is the normalization 
volume of the phonon modes, and cr e/ / h are the deforma- 
tion potential constants for electrons and holes. 

The Hamiltonian H sc [Eq. has the structure of 
the independent boson model and may be diagonalizcd 
exactlySS. To this end, one defines the operator 



W = exp 



^2ala n gl(k)b k - H.c. 



where g n {k) — f n (k)/ (hujk), and writes the Hamiltonian 
Hq in terms of the operators 



a n = W^aJW = a n W n 



(4) 



and 



where 



f3k=W% k W = b k + J2 a i a ng n (k), (5) 



W n = exp 



5>;(fc)&fc-H.c. 



In the calculations one uses the fact that [W n , W m ] 
for non-overlapping exciton states. The result is 
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where E n = e n — J^ fc \f n (k)\ 2 / (fnok)- The new operators 
(@J and © represent excitons in the dots surrounded 
by coherent lattice displacement field and phonon modes 
shifted by the presence of a charge distribution. 

We will associate the value |0) of the logical qubit with 
the state of the QD array in which a certain subset oiN/2 
dots, labelled by indices m, . . . un/2, is occupied by exci- 
tons and the value |1) with the state in which the other 
N/2 (n 1; . . . n' N/2 , with n ; ^ n' v for allZ, V = 1, . . . , N/2) 
dots are occupied. Thus, the basis states of the logical 
qubit are 
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where |g) is the ground state of the system (with all dots 
empty). In our discussion we will assume that a control 
field is available that couples the relevant pair of states of 
the multi-QD register and induces a transition without 
leaving the subspace spanned by these states. In partic- 
ular, we assume the availability of a control Hamiltonian 
of the form 



He = \f{t)e 



-iAEt/h f 



H.c. 



(7) 
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which transfers N/2 excitons from one subset of dots to 
the other. Here fit) is the envelope of the control field 
and 



N/2 

AE = ^2E n 

1=1 



N/2 

1=1 



is the difference between the total energies of the logical 
states |0) and A practical implementation of such 
a coupling is a separate issue. For transferring a single 
exciton between two dots, a realization of an effective 
Hamiltonian of this kind has been proposed, using the 
Forster coupling and the optical Stark eSecijS^^. 

In terms of the redefined degrees of freedom [Eqs. (0} 
and 10] , the control Hamiltonian @ to the leading order 
in phonon coupling reads 
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with 



G(k) = [g n {(k)+...+g n ' N/2 {k)}-[9nAk)+---+9n N/2 (k)]. 

The phonon perturbation in Eq. (JSJ leads to reservoir 
excitations accompanying any operation on the carrier 
subsystem but does not drive this subsystem outside the 
subspace spanned by the logical qubit states. 

In our calculations the wave functions of confined elec- 
trons and holes will be modelled by identical Gaussians, 
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where the parameters I, l z are related to the QD size in 
the xy plane and along the z direction, respectively. For 
simplicity, we assume here that the wave functions in all 
QDs are identical. As shown in the Appendix, including a 
small variation of dot sizes and relative electron and hole 
confinement widths will only lead to inessential quanti- 
tative corrections. For this choice of wave functions, the 
formfactor [Eq. (J3J] has the explicit form 
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where k± = (k% + k 2 ) 1 / 2 . 

We will assume that the operation is performed by a 
Gaussian control pulse, 
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where a is the angle determining the gate (corresponding 
to the angle of rotation on the Bloch sphere) and t p is 
the gate duration. 

In Tab.^the material parameters (corresponding to an 
InAs/GaAs system) are collected. 



Deformation potential coupling a c — <Jh 8 eV 



Crystal density p c 
Speed of sound (longitudinal) c 
Long-time decoherence times- to 
at T = K 
at T = 10 K 



5360 kg/m 3 
5150 m/s 

2 ns 
0.5 ns 



TABLE I: System parameters used in the calculations. 



III. PHONON PERTURBATION DURING A 
QUANTUM GATE 

In the interaction picture with respect to H sc the sys- 
tem evolution is generated by the Hamiltonian 



H = \}{t)a x + ~f(t)a y Rit), 



(11) 



where 



R(t) =i^G(fe)(e" 
k 



and a x<y are Pauli matrices. 
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The evolution in the absence of phonon perturbation 
is described by the unitary operator t/o(i) generated by 
the first term in Eq. I|ll|) . 
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The effect of the phonon coupling on the driven dy- 
namics may be calculated using the second-order (Born) 
expansion of the evolution equation for the density 
matrijiiS'Si, including the first term in Eq. (|llfl exactly 
and the second one as a perturbation. The reduced den- 
sity matrix of the qubit is written as 



p{t) = U Q it)[p Q +p^it)]Ulit), 



(12) 



where the correction it) is calculated from a pertur- 
bation expansion 
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Here Ttr is the trace with respect to the reservoir 
(phonon) degrees of freedom and V(t) — Uq (t)Vf7o(t), 
where V is the second term in Eq. (|11|) . We will as- 
sume the initial state of the compound system in the 
form go = Po ® Pt, where pt is the thermal equilibrium 
state of the phonon bath and the initial state of the qubit 
subsystem is pure, po = l^'oXV'ol- 
Defining the operator 
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and the phonon spectral function 

R(lu) = -L J dt(R(t)R)e lu \ 
and representing the Hcavisidc function as 



e(t) 

one may write 

-,( 2 
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^(t) = -i [h t ,p ] - - {A t , Po } + <f> t [p ], (15) 
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FIG. 1: The spectral characteristics of the driving, S(ui), for 
two rotation angles on the qubit Bloch sphere and for a series 
of pulse durations as shown, for Gaussian pulses [Eq. llUH l. 
The scales on the two figures are different. 
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with 

A t = J dwR(u)Y^(u)Y(w), 

f , „, n« f du' YUu')Y(uj') 

h t = / duR(u>) V j V • 

J J 2ir ui' — oj 

where V denotes the Cauchy principal value. 
The explicit form of the operator Y(ui) is 



Y(u)= t -F(u,)\+)(-\+ i -F*(-u,)\-){+\, 

where |±) = (|0)± |1))/V2 and F(u) is a nonlinear spec- 
tral characteristics of the control field, 

1 r°° 

F(u) = - / dte^ftye**®. (19) 

™ J -co 

We assumed here that the initial time to is before the 
control field has been switched on and the final time t is 
after it has been switched off so that the limits of inte- 
gration in the above formulas could be extended to ±oo. 

To measure the effect of the phonon-induced dephas- 
ing on a quantum gate we will use the error S, defined as 

S — 1 ~ F 2 , where F = (V>o|^ t (co)p(oo) L / (co)|V>o) 1/2 is 
the fidelity distance measure^ between the ideal (unper- 
turbed) final state ^(oo^oX^ol^o (°°) and the actual 
state p(oo). Substituting Eq. (|12fl into this definition one 
finds 

<5 = -^o|p (2) (oo)|^o)= / dwJJ(«)|«|y(w)|^)| 2 , 



where we used the explicit expression (|15f) and \ipQ~) is a 
state orthogonal to j^o) in the logical qubit space. Choos- 
ing 



4>o = sin-|+> -e^cos-|-), 
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one gets 

|(V^|rMhM 2 = 



F(u)e iv cos 2 - - F*(~uj)e-' lip sin 2 - 



Since, in general, the qubit may initially be in any state, 
we average the error over the Bloch sphere of initial 
states, i.e., over the angles (9,<p), which yields 

S(co) = \(^\Y(^ )\l = ±(\F(u,)\ 2 + |F(-c)| 2 ). 

Since this is even in u, the error can be written in the 
form 

,S= / ( /., ( -oth| J^L] ^-S(u>), (20) 



where 



J(u) 
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is the standard spectral density of the phonon reservoir. 

From Eq. (|19|l it can be seen that F(u>), and therefore 
also S(u>), show simple scaling behavior with respect to 
the duration of the control pulse. Namely, for a family of 
pulses of a fixed shape and area, f(t) = (Ha / tq) f (t / tq) , 
where f(x) is the fixed shape of the pulse and To is the 
pulse duration, the spectral function scales as S(uj) = 
S(toiv) (see Fig. P). Therefore, for sufficiently long con- 
trol pulses the error is determined by the low frequency 
part of the phonon spectral density. To be specific, it 
can be seen from Eq. 1121)1) that if, for low frequencies, 



J(u>) ~ LJ n then 5 ~ r t 



-ra+l 



or S ' 



-n+2 



for long enough 



pulses at low and high temperatures, respectively^. Ob- 
viously, slowing down the driving in order to reduce the 
error is possible only to a certain extent, since for longer 
durations the error becomes dominated by decoherence 
effects related, e.g., to finite exciton lifetime and phonon- 
assisted transitions, which accumulate with time (linearly 
for short times). The interplay of these two sources of 
dephasing creates a trade-off situation and sets a lower 
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FIG. 2: The two layouts of the QD array: 3-dimensional 
(left) and linear (right). The dark dots are occupied in the 
logical |1) state and the white dots in the |0) state. In the 
linear layout, the logical qubit with 2N dots is constructed 
by interchanging occupied and empty dots in the N dot qubit 
and appending it to the original one. 



bound on the total error achievable for various pulse 
durations 17 . 

Although, in general, the phonon spectral density for a 
single exciton confined in a QD depends on the dot geom- 
etry, the low frequency behavior is only material-specific 
(see Appendix) and can only be modified by engineer- 
ing the properties of phonon mode a 20 i 21 . However, as we 
show in the next Section, the low-frequency behavior of 
the effective phonon spectral density may be modified by 
collective encoding in an array of QDs, which allows one 
to reduce the total error. 



IV. ERROR REDUCTION BY COLLECTIVE 
ENCODING 

We consider two kinds of arrays of quantum dots. The 
first one, referred to as 3-dimensional (or 3D) , is formed 
by stacking 2 QDs along the strong confinement (z) axis, 
arranging 4 QDs in a rectangle in the xz plane or placing 
8 QDs in the corners of a cuboid (Fig. [21 left). In the sec- 
ond layout, which we will refer to as linear, N = 2 n dots 
are arranged along the z axis, as shown and explained in 
the right part of Fig. [2J 

For N = 2" quantum dots in the 3D layout one finds 
the coupling constants G n (for an appropriate choice of 
origin) 



k z d z 



Gi(fc) = 2 lS in^G (fc), 
G 2 (fc) = (2z) 2 sin^sin^G (fc), 
G 3 (fc) = (2z) 3 sin^sin^sin^G (fc), 
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FIG. 3: The phonon spectral densities for the physical ex- 
citonic qubit (solid line) and for logical qubits composed of 
N physical qubits (dashed lines) for two kinds of QD arrays: 
3-dimensional (a,b) and linear (c,d) with geometrical param- 
eters as indicated. 



is the coupling constant for a single dot placed at the 
origin and the distance parameters d, d z are defined in 
Fig. [21 From Eq. H21fl one then finds the phonon spectral 
densities for the logical qubits in the low frequency sector, 
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These spectral densities are plotted in Fig.[^a,b). 
For the linear array of QDs one finds 

n-2 

G n (fc) = (2i) n J] sin(2^ z 4)G (fe), (23) 

i=-i 



which leads to the spectral density J\(w) as above and 
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shown in Fig. E£c,d). 

In both layouts, J n (u>) ~ ui 2n+3 for small u so that for 
larger logical qubits the spectral density vanishes more 
quickly for to — > 0, which will result in a faster decrease 
of the pure dephasing error with growing duration of the 



6 



10" 



10" 



10" J 
10" 4 

10° 
p 10" 1 
10 2 
10" 3 







a = nil 

i \ \ \ 
(a) Is \ 


a - Till 

N=l — % 

N-2 — \ 

" N=4 v\ 

(b) N=8 - ^\ 


0.1 1 0.1 1 


/Z\ 7-0 K 
/.•;. — a = 5nll 

~^ V\ j 
(0 , % \ 


,v-N / - 10 K 
. , a - 5.T/2 

(d) , . \\\. 



0.1 



10 



0.1 



?p [ps] 



10 



FIG. 4: The phonon-induced dynamical error for a 7r/2 and 
5tv/2 rotation on a physical qubit (solid line) and on a 3D 
logical qubit with 2,4, and 8 dots (dashed and dotted lines) 
with sizes I — 4 nm, l z = 1 nm and separations d — 15 nm, 
d z — 4 nm, at two different temperatures as shown. 



control pulse, as discussed in the previous Section. In 
order to see this, we calculate the phonon-induced er- 
ror using Eq. (|20(l with the appropriate spectral density 
J n (w) [obtained from Eq. H21|l , for arbitrary uj] for vari- 
ous logical qubits. 

Let us start the discussion with a logical qubit in the 
3D layout, composed of QDs of dimensions typical for 
self- assembled systems (Fig. 0J. First, let us consider 
the 7r/2 rotation of the qubit (i.e., the single-qubit gate 
e -ma y /4y p or m trafast pulses (t p < 100 fs), in most 
cases the error grows with the number of QDs in the 
array. This regime corresponds to very broad spectral 
functions S(ui) (see Fig. [IJ, so that the whole range of 
phonon modes contribute to the dephasing. Since con- 
siderable destructive interference of phonon excitations is 
unlikely to appear over the whole wide frequency sector 
of the phonon spectral density, it may be expected that 
driving a number of dots simultaneously should lead to 
more excitation of the phonon reservoir and to more de- 
phasing. In terms of the spectral densities (Fig. 01, this 
effect is reflected by the growing overall magnitude of 
J(u). 

This situation changes for longer pulses (r p ~ 1 ps). 
Now only the long wavelength phonons contribute to de- 
phasing. From the general discussion presented above 
it is clear that for sufficiently long pulses large collec- 
tive qubits will lead to lower phonon-induced error than 
smaller ones due to more favorable scaling of the error 
with the pulse duration. Indeed, in Fig. 0] the value of 5 
decreases for long enough pulses (r p > 1 ps) according to 
a power-law dependence with higher negative exponents 
for larger arrays. 

It should be stressed that the reduction of dephasing 
achieved with collective encoding is not just due to ex- 
tending the effective system size and the resulting weak- 
ening of the overall phonon couplingii^. Such an ef- 
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FIG. 5: Comparison of the dynamical error for three different 
geometries of the QD array: small, closely spaced dots in the 
3D layout (a,b), flat dots of moderate size (like in Fig. [IJ in 
the linear layout (c,d) and large dots in the linear layout (e,f). 



feet might quantitatively reduce dephasing but could not 
change the exponent of the low-frequency behavior of 
J(u>). Moreover, the reduction of J(u>) at u> — > appears 
only for the special definitions of logical qubits shown in 
Fig. [3 For a different assignment of physical QDs to the 
logical qubit states, the dephasing may even be increased, 
which shows that the effect is actually due to destructive 
interference between the reservoir excitations caused by 
the simultaneous switching of the appropriately arranged 
physical qubits. 

For large angles a, corresponding to full Rabi rota- 
tions around the Bloch sphere, the error dependence on 
the pulse duration becomes more complicated due to the 
presence of many maxima in S(ui) (Fig. [JJ. While the 
power-law decrease of 6 at r p > 1 ps corresponds to the 
situation when all of S(u>) lies in the low-frequency sector 
of J(u>), the minimum at r p ~ 100 fs appears when the 
area of large J(ui) coincides with the local minimum be- 
tween the maxima of S(u). Since the collective encoding 
improves only the low-frequency behavior of J(lu) only 
the former is advantageously affected by increasing the 
size of the QD array. 

In the ultrafast limit, the dephasing of a single physi- 
cal excitonic qubit increases for decreasing dot sizes, since 
the 1/2 momentum cut-off in the formfactor [Eq. (JSJ] is 
shifted towards larger momenta and the total area of the 
phonon spectral density grows. This can also be seen in 
Fig- Gib) > corresponding to QDs with the lateral size re- 
duced to I = 1 nm. For longer pulses, this difference in 
J(u>) becomes less important since only the low-frequency 
part of this spectral density is relevant. Indeed, com- 
paring Fig. |SJa,b) (corresponding to smaller dots) with 
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Fig.|3Ja,b) one notes that the error for ultrashort pulses is 
much larger for smaller dots but in the region of r p > 1 ps 
the error for a single dot (solid line) grows very little when 
the dot size is reduced. On the other hand, smaller dots 
allow smaller distances d between the physical qubits in 
the array. As a result, the modulation of J(lu) resulting 
from the multi-QD interference attains a longer period 
Aw ~ c/d and the sector of low frequencies, where the 
spectral densities J„(w) are ordered according to their 
asymptotic power-law behavior, extends [see Fig. E[b)]. 
Now, the improved behavior of the spectral densities is 
fully reflected by the reduced values of errors for growing 
number of QDs in the array already for r p > 300 fs. 

The linear layout allows the dots to be tightly stacked 
along the direction of smallest size (growth direction of 
the semiconductor structure) but, on the other hand, re- 
sults in a growing distance between the physical qubits on 
the opposite ends of the chain. This leads to a quickly 
growing pre-factor at the low-frequency power-law ex- 
pression [compare Eqs. (|22a|l and l|24ajl ]. As a result, the 
spectral densities for larger arrays cross Ja(w) at very low 
frequencies (see insets in Fig. [2J and one needs longer 
pulses to take advantage of the lower values of J n (co>) for 
n > 1, as can be seen in Fig. Otc,d). 

Obviously, the quantitative results presented here de- 
pend on the interplay of the physical coupling parameters 
and the wave number dependence of the original (physi- 
cal) coupling in the long wave length limit. In any case, 
however, the spacing between the dots is crucial. This is 
illustrated by Fig.|3Je,f), where the total error for a linear 
array of large and widely spaced dots is plotted. From 
the spectral density shown in Fig. EJd) it is clear that 
such a geometry leads not only to the overall narrowing 
of the phonon spectral characteristics (due to large size) 
but also to the squeezing of the region when the functions 
J(u>) show the advantageous low-frequency behavior. As 
a result, increasing array size leads to error reduction 
only for very long pulses. 

In a real system, the dynamically induced pure de- 
phasing is not the only error source. Additional error 
comes from the finite exciton lifetime and other pro- 
cesses, including phonon-assisted transitions to excited 
exciton statesi2iM, dephasing by phonon scatteringiSSi 
and probably other effects whose exact nature seems to 
be poorly understood. The common feature of all these 
dephasing sources is that they can be described within 
a Markovian model leading to an exponential decay of 
coherence with a single time constant tq which is ex- 
perimentally accessible via exciton line broadening. For 
r p <C To their contribution to the overall error grows lin- 
early with the gate duration t p . Therefore, we add a 
contribution of Sq = t p /tq to the total error. For brevity, 
this additional contribution will be referred to as related 
to "finite lifetime". Depending on the physical process 
underlying this error, the time constant To may depend 
in various ways on the number of physical qubits in the 
array. Here, we will consider two scenarios. In the first 
one, To is independent of N (which assumes some kind 
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FIG. 6: The total error, including a contribution from finite 
exciton lifetime and other long-time scale decoherence pro- 
cesses for typical flat dots (with sizes as shown) in 3D (a,b) 
and linear (c,d) layout at two different temperatures. 



of collective action of the dephasing), while in the sec- 
ond one it is proportional to 1 /N (which corresponds to 
individual decoherence of the physical qubits). 

Reaching the low frequency regime, necessary for re- 
ducing the dynamical error, requires long control pulses 
and leads to growing values of <5o- As a result, for long 
pulse durations the error is always dominated by the 
linearly growing finite lifetime contribution. The inter- 
play of the dynamical phonon-induced error (growing for 
short control pulses) and the finite lifetime decoherence 
(increasing for long pulses) leads to well-defined driving 
conditions where the error has its minimumii. As shown 
in Fig. El in the first model of decoherence (with lifetimes 
To = 2 ns and 0.5 ns at T — K and 10 K, respectively^,) 
the optimal total error is reduced by a factor of 2 at T = 
K and almost by a factor of 4 at T = 10 K. It is interest- 
ing to note that for these life times large logical qubits 
N > 2 bring very little improvement over the simplest 
one (N = 2) and only in the 3D layout. 

For the second decoherence model, the reduction of 
the effective lifetime of the array considerably restricts 
the possible error reduction. At T — one obtains only 
a fidelity gain of 17%. More error reduction (by a factor 
of 2) is achievable at higher temperatures. Due to large 
increase of the error with growing number of physical 
qubits, only the encoding in N = 2 physical qubits brings 
the desired effect. 

Thus, nanosecond decoherence times observed in the 
currently manufactured structures strongly restrict the 
degree of reduction of dephasing achievable with collec- 
tive encoding, especially in the independent decoherence 
model. For realistic sizes and separations of the QDs in 
the array, taking advantage of the reduced spectral den- 
sity requires pulses of durations exceeding 1 ps, which 
immediately leads to an error of order of 10 -3 , further 
growing for increasing number of dots in the array. How- 
ever, rapid progress of QD engineering may lead to sys- 
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FIG. 7: As in Fig.||j]but with decay probability proportional 
to the number of physical qubits. 
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FIG. 8: The total error for a 3D logical qubit with geometry 
as in Fig. [Tfa.bl for (hypothetical) long characteristic times 
of the exciton lifetime, with decay-related error proportional 
to the number of physical qubits, for a = n/2 at two temper- 
atures. 



tems with considerably longer coherence times, as sug- 
gested by the observation of long-lived excitons in par- 
ticular QD systems 37 . An advantage of the collective 
qubit is that such reduction of the long time decoherence 
is much more effectively reflected by the reduction of the 
total error. As can be seen in Fig. extending the ex- 
citon lifetime leads to a reduced error both for a single 
physical qubit and for a collective qubit. However, the 
error reduction is much larger for the logical qubit. For 
instance, the 2-QD encoding (Fig. Q) at T — reduces 
the error by 17% for To = 2 ns, which grows to 31% if the 
lifetime is extended to 5 ns and to a factor of almost 1.7 
for To = 30 ns. There is even more gain at higher temper- 
atures. For T = 10 K, one obtains an error reduction by 
a factor of 2.8 if the lifetime is extended to 2 ns (which is 
the actual radiatively limited value in some structures 9 ) 
and by almost an order of magnitude for tq = 30 ns. 



CONCLUSIONS AND OUTLOOK 



The discussion presented above shows that engineering 
of the effective phonon spectral density by collective en- 
coding of a logical qubit in an array of physical (excitonic) 
qubits may lead to a reduction of the phonon-induced de- 
phasing during control operations on the qubit. This is 
possible since for long enough pulses (slow driving) only 
long- wavelength part of the phonon reservoir is involved 
in the dephasing dynamics which extends the effective 
correlation length of the reservoir and opens the possi- 
bility for collective encoding. Such collective encoding 
reduces the total error resulting from the joint action 
of various dephasing mechanisms and increases the gain 
from possible further extension of the exciton life times. 
With the collective encoding it seems feasible to reduce 
the total error for a single-qubit operation to values of 
order of 10 -3 even at the relatively high temperature of 
T = 10 K, which is remarkable in view of the recent 
progress in the quantum error correction theory under 
realistically strong dephasing 3 ^. 

The present analysis was restricted to Gaussian pulses. 
It is known, however, that pulse shaping^ or pulsed con- 
trol techniques^ lead to reduced phonon dephasing ef- 
fect. The reach structure of the effective phonon spec- 
tral density for collective qubits suggests that combining 
the pulse shaping techniques with the collective encoding 
may be very fruitful. Another possibility of reducing the 
overall error might be to combine the collective encoding 
against dynamical dephasing with the noiseless (subradi- 
ant) encoding against radiative decay^. 

In more general terms, our discussion shows that cer- 
tain details of the phonon spectral density may be of es- 
sential importance for the dephasing of a system driven 
by finite pulses. Obviously, the system response to an 
ultrashort pulse, as observed in optical experiments 10 , 
yields valuable information on the overall strength of the 
phonon-induced dephasing. However, we have seen that 
the degree of decoherence for finite pulses is governed 
only by the low-frequency part of the relevant spectral 
function and may decrease even though the dephasing in 
the ultrafast limit increases (e.g., for arrays of smaller 
but more closely spaced QDs). It is therefore possible 
that systems based on various materials 3 ^ and restricted 
geometries^ - may offer previously unnoticed possibilities 
in the slow driving regime. Understanding of the pecu- 
liarities of the carrier-phonon dynamics underlying the 
dephasing of driven systems may also be helpful for de- 
signing robust schemes for the optical control of spin 
qubits. 
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port. 
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APPENDIX A: UNIVERSALITY OF 
CARRIER-PHONON COUPLING IN THE LONG 
WAVE LENGTH LIMIT 



The formfactor can now be written as 



In this Appendix we show that for small variation 
of individual dot sizes and for small difference between 
the electron and hole localization in a dot, the exciton- 
phonon coupling constant can be written in the form J2J, 
with the formfactor replaced by Eq. © which depends 
only on the average width of the electron and hole wave 
functions in the QD array. 

In general, the electron and hole wave functions in the 
nth dot may be characterized by different localization 
widths l„ (for simplicity, we still assume that the ex- 
citon wave function is a product of Gaussians isotropic 
in the xy plane and with a fixed width l z ). The carrier- 
phonon interaction is then described by Eq. J2J withSi 

T n {k) = a oe - fc i^ >2 /4-^/4 _ z he -kH^/4-klll/4^ 

where CT G /h = o' c /h/( cr c — &h)- We define the averaged 
wave function width 



i» = JL 

27V 



£(4 e) ) 2 + BW 



and the deviations from the average 



(A4 c/h) ) 2 = (e /h) ) : 



-I 2 



x <r e e 



kl(Al^f/4_~ he -kl(Al^f/4 



It is clear that kj_ is cut off at the value of ~ l/l, hence 

k±Aln^ h ^ 1 if Al„ *C I, i.e., when the variation 
of dot sizes is relatively small. The exponents may be 
therefore expanded into power series, 



T n (k) 



_ -k 2 ± l 2 /4-k 2 J 2 z /4: 



(e)N2 



1-5, 



k 2 ± (Al 



(h)\2 



where we took into account that a c — Ch = 1- Retaining 
only the leading term one arrives at Eq. ©. It should 
be noted that independently of the shape of the wave 
functions, normalization condition requires that for long 
wavelengths T n (k) = 1 + 0(fc 2 ), which is sufficient for 
deriving Eqs. (|22at c) and I24al b). 
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